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NOMENCLATURE 


Q  2nd  order  tensor  (shear)  of  velocity  gradient,  s'2 
rc  vortex  core  radius  in  the  PIV  image,  m 
R  blade  radius,  m 

U  local  flow  velocity  in  the  horizontal  direction,  m/sec 

V  local  flow  velocity  in  the  vertical  direction,  m/sec 

W  local  flow  velocity  in  the  out-of-plane  direction,  m/sec 

as  shaft  angle,  positive  aft,  degrees 

Tv  vortex  circulation,  m2/s 

yz  vorticity  in  the  out-of-plane  direction,  s'1 

X2  2nd  order  tensor  (strain)  of  velocity  gradient,  s'2 

§  elastic  torsion,  degrees 

Q  rotor  operating  speed,  radians/sec 

<nz  angular  velocity  (=  yz/2),  s'1 

ax?  y  standard  deviations  of  the  vortex  centers  in  the  x-  or  y-direction 
0X  PIV  image  rotation  in  the  horizontal  axis,  degrees 

0y  PIV  image  rotation  in  the  vertical  axis,  degrees 

\j/  rotor  azimuth,  degrees 

ACRONYMS 

AFDD  Aeroflightdynamics  Directorate 

BL  Baseline 

BVI  Blade  Vortex  Interaction 

CA  Conditional  Averaging 

DLR  Deutsches  Zentrum  fur  Luft-und  Raumfahrt,  German  Aerospace  Center 

DNW  German-Dutch  Wind  tunnel 

Dpt  Data  Point 

LLF  Large  Low-speed  Facility 

HART  Higher-harmonic  control  Aeroacoustic  Rotor  Test 

NASA  National  Aeronautics  and  Space  Administration 

ONERA  Office  National  d'Etudes  et  de  Recherches  Aerospatiales 

PIV  Particle  Image  Velocimetry 

Pos  Position 

SA  Simple  Averaging 

3-C  3 -Component 


ABSTRACT 

Since  the  completion  of  the  DNW  test,  HART  II  partners  have  put  forth  significant  efforts  to  derive 
vortex  properties  from  the  measured  data.  These  properties  include  vortex  center,  rotation  angles  of  the 
vortex  axis  in  the  PIV  measurement  plane,  maximum  vorticity,  core  radius,  and  swirl  velocity. 
Differences  among  the  partners6  results  were  found  for  some  of  these  vortex  properties,  which  required 
more  efforts  to  assess  the  physical  significance  of  the  differences.  This  was  because  the  derived  vortex 
properties  were  an  integral  of  all  the  sub-processes  in  the  complex  PIV  data  analysis  procedure,  and  more 
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information  was  needed  to  resolve  the  differences  in  vortex  properties.  The  objectives  of  this  report  are  to 
discuss  and  summarize  the  parameters  and  techniques  that  could  contribute  to  these  differences.  The 
detailed  comparisons  in  this  report  have  helped  establish  a  strong  foundation  for  a  PIV  data  analysis 
methodology  in  determining  vortex  properties. 


INTRODUCTION 

Particle  Image  Velocimetry  (PIV)  measurement  data  can  provide  a  significant  amount  of 
information  about  vortex  core  dynamics  in  a  complex  flow  field.  There  have  been  many  reports  about  the 
PIV  data  analysis  techniques  applied  to  the  rotor  vortex  core  properties  [1-7],  but  the  techniques  have  not 
been  standardized  in  the  rotorcraft  field,  which  was  partially  due  to  the  complexity  of  the  inherent, 
rotational  flow  characteristics.  Rotor  wakes  interact  with  rotor  blades  and  occasionally  with  returning 
vortex  wakes,  which  is  called  Blade -Vortex  Interaction  (BVI).  This  BVI  phenomenon  can  be  typically 
seen  in  descending  flight. 

One  of  the  difficulties  in  deriving  vortex  properties  from  3-C  PIV  measurements  is  that  the  plane 
normal  to  the  vortex  axis  is  unknown  a  priori.  Another  difficulty  is  that  the  convective  velocity  of  a 
vortex  is  also  unknown  since  it  is  a  sum  of  a  free  stream  velocity  and  the  velocity  induced  by  other 
vortices.  In  addition,  a  vortex  is  not  always  rotationally  symmetric;  often  a  sheet  of  trailed  vortices  is 
found  spiraling  around  a  vortex. 

This  report  describes  the  efforts  of  an  internationally -joint  HART  II  team  for  extensive  analyses  of 
the  PIV  database  that  was  acquired  during  a  wind  tunnel  test  of  the  HART  II  rotor.  A  brief  description  of 
the  PIV  database  is  also  provided,  and  the  techniques  for  deriving  vortex  properties  such  as  vortex  center, 
rotation  angle  of  the  vortex  axis,  maximum  vorticity,  core  radius,  and  swirl  velocity  are  followed. 
Differences  in  resulting  vortex  properties  among  the  partners  are  addressed  and  discussed.  Specific 
findings  from  this  study  are  provided  in  the  conclusions. 


HART  II  3-C  PIV  DATA 

In  2001,  an  international  cooperative  test,  called  the  Higher-Harmonic-Control  Aeroacoustic  Rotor 
Test  (HART  II)  was  conducted  by  the  US  Army  Aeroflightdynamics  Directorate  (AFDD),  Deutsches 
Zentrum  fur  Luft-und  Raumfahrt  (German  Aerospace  Center,  DLR),  Duits-Nederlandse  Windtunnel 
(German-Dutch  Wind  Tunnel,  DNW),  National  Aeronautics  and  Space  Administration  (NASA),  and  the 
Office  National  d‘ Etudes  et  de  Recherches  Aerospatiales  (Onera)  at  the  large  low-speed  facility  (LLF)  of 
the  DNW  wind  tunnel.  Primary  objectives  of  the  test  were  to  measure  rotor  wake  information  extensively 
using  3 -component  (3-C)  PIV,  blade  surface  pressure  distributions,  blade  deformations,  and  acoustic 
signatures  [8-10]. 

The  rotor  was  a  40%  Mach-scaled,  hingeless  Bo  105  main  rotor  (HART  II),  and  had  a  precone  of 
2.5  degrees  at  the  hub,  operating  at  1041  rpm.  The  blade  was  2  meters  long  with  -8  degrees  twist,  and 
consisted  of  a  NACA-23012  airfoil  section  with  a  trailing  edge  tab.  The  chord  length  was  0.121  meters. 
The  hover  tip  speed  was  218  m/s  with  a  tip  Mach  number  of  0.64.  The  rotor  shaft  angle  in  the  test  was  5.3 
degrees  aft. 

Figure  1  shows  the  HART  II  3-C  PIV  setup  that  included  three  double-pulsed  Nd:YAG  lasers  and 
four  digital  cameras  mounted  onto  a  remote-controlled  traverse  system.  This  traverse  system  was  located 
outside  of  the  shear  layer  in  the  open  jet  configuration  and  moved  in  the  direction  of  flow.  Four  cameras 
on  the  traverse  system  were  used  to  collect  3-C  PIV  data.  One  pair  of  cameras  used  100mm  lenses  for  a 
larger  observation  area  (DNW  measurement  window)  and  the  other  camera  pair  used  300mm  lenses  to 
investigate  the  tip  vortex  with  high  spatial  resolution  (DLR  measurement  window). 

The  PIV  measurement  plane  coordinate  system  has  the  positive  x-axis  to  the  right  direction,  the 
positive  y-axis  up  and  the  positive  z-axis  out  of  plane.  The  rotor  hub  coordinate  system  is  aligned  with  the 
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rotor  shaft,  such  that  the  positive  x  is  toward  the  tail  boom,  the  positive  y  is  toward  the  right  wing,  and  the 
positive  z  is  upward  along  the  shaft.  The  wind  tunnel  coordinate  system  is  defined  as  the  inertial  frame 
without  considering  a  shaft  tilt:  the  positive  x  upstream,  the  positive  y  to  the  right  (pilot1 2 3 4 5 6 s  view),  and  the 
positive  z  down. 

Figure  2  shows  the  location  and  orientation  of  the  PIV  measurement  planes  for  the  HART  II  test. 
To  avoid  interference  of  the  blades  with  the  measurement  plane,  the  reference  blade  was  located 
alternatively  at  an  azimuth  of  20°  or  70°.  In  each  measurement  plane  (black  line  segment),  one  hundred 
individual  PIV  images  were  obtained,  and  the  numerals  represent  the  identification  number  of  each  PIV 
measurement  position.  The  measurement  planes  were  oriented  at  an  angle  of  149.35°  for  the  advancing 
side  and  30.06°  for  the  retreating  side  with  respect  to  the  positive  hub  x-axis. 

Figure  3  shows  an  individual  image  for  Dpt  995  (position  21),  taken  by  the  camera.  The  upper 
image  (0.45m  x  0.38m)  is  the  global  view  taken  from  the  DNW  camera  with  the  apparent  center  of  swirl 
motion  (vortex  center)  in  the  flow  field.  The  close-up  view  image  (0.15m  x  0.13m)  from  the  DLR  camera 
shows  the  indicated  area  of  the  vortex  with  a  higher  resolution.  However,  the  high  resolution  image  does 
not  indicate  a  clear  swirl  motion. 

The  velocities  (U,  V  and  W)  and  the  vorticity  for  image  #1  at  the  position  17  are  shown  in  Fig.  4. 
The  vorticity,  yz  was  derived  using  the  relation: 


yz  =2  coz  =  k  (curl  V)  = 


dV 

dx 


8U 

dy 


(1) 


where  ooz  is  the  angular  velocity  and  V  =  (U,  V,  W).  In  the  figure,  the  peak  vorticity  and  noise  of  the 
measurement  are  easily  observed.  The  noise  may  result  from  natural  turbulence  and  unsteadiness  of  the 
flow  field,  and  the  vortex  center  could  be  affected  by  rotor  test  stand  vibration,  camera  support  vibration 
and  other  unknown  resources. 


VORTEX  CENTER  AND  CONDITIONAL  AVERAGING 

Averaging  a  large  number  of  velocity  vector  fields  will  reduce  the  noise  inherent  in  an  individual 
image.  It  is  straightforward  to  overlay  all  of  the  individual  velocity  fields  in  the  PIV  images  and  then  to 
average  them  to  a  single  velocity  field.  This  is  called  simple  averaging  (SA)  technique  and  works  well 
when  the  vortex  centers  are  mostly  collocated  from  one  image  to  another.  In  the  test  setup,  the  rotor 
model  stand  was  attached  to  the  tunnel  sting.  The  sting  and  the  camera  supports  were  slightly  vibrating  at 
low  frequencies  due  to  elasticity.  The  test  rig  vibrations  were  transferred  to  the  blades,  which  resulted  in 
scattering  of  the  vortex  centers  between  individual  PIV  images.  In  addition,  natural  instabilities  of  helical 
vortex  systems  added  vortex  wandering  with  the  vortex  age.  This  situation  makes  it  necessary  to  conduct 
a  conditional  averaging  (CA)  technique  in  order  to  prevent  smearing  of  the  vortex  which  would  occur  if 
the  SA  technique  is  employed.  The  details  of  a  conditional  averaging  technique  are  given  in  Ref.  [11]. 
Figure  5  shows  the  procedure  for  conditional  averaging  of  velocity  fields,  which  is  briefly  given  as 
follows: 

1)  Specify  the  sense  of  rotation  of  the  vortex  under  investigation,  since  the  vortices  in  the  opposite 
rotation  could  be  in  the  same  image. 

2)  Identify  the  vortex  center  location  in  each  PIV  image  using  a  sub-grid. 

3)  Compute  the  standard  deviation  of  vortex  centers  for  all  100  PIV  images. 

4)  Within  a  user-specified  Z  score  (Z  ~  1.5)  select  80-95  individual  images  (geometric  outlier 
removal).  Re-compute  the  mean  vortex  center  position  based  on  selected  PIV  images.  Note  that  Z 
score  is  the  sample  score  non-dimensionalized  with  a  standard  deviation. 

5)  Perform  a  conditional  averaging  by  collocating  selected  images  with  individually  computed 

vortex  centers.  The  vortex  center  of  each  image  does  not  necessarily  lie  on  a  grid  point,  which 
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requires  re -meshing  of  the  individual  images  and  then  interpolating  their  velocities  in  the  re¬ 
meshed  images. 

A  vortex  can  be  defined  using  a  various  form  of  velocity  tensors.  We  consider  here  four  different 
velocity  gradient  tensors  for  identification  of  a  vortex:  vorticity  yz,  Q  criterion,  X2  and  convolution  of  X2 
[11-12].  The  first  order  tensor,  vorticity,  is  widely  used  in  identifying  a  vortex.  Higher  order  tensors  are 
understood  to  display  a  clearer  image  of  a  vortex  than  the  lower  order  tensor.  The  second  order  tensors 
include  Q  criterion,  X2  and  convolution  of  X2.  The  first  two  definitions  are  similar  to  each  other,  and  their 
explicit  expressions  for  the  two-dimensional  velocity  field  are  given  as  follows  [11,  13-14]: 


"et/f 

ydx  ) 


rdV^2 


dy 


dV  8U 
dx  dy 


'du  avf 

v  dx  +  dy  , 


dU_dV_  dV  dU 
dx  dy  dx  dy 


]_ 
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dU  dV  |  dV  dU 
dx  dy  dx  dy 


(2) 


Note  that  vorticity  (the  first  order  tensor,  yz)  provides  information  on  the  sense  of  rotation  via  its  sign,  but 
the  second  order  tensors  (Q  and  X2)  do  not  provide  this  information  —  the  sign  of  their  peak  values  is 
independent  of  the  sense  of  rotation.  The  convolution  of  X2  maps  the  X2  distribution  into  a  bell-shaped 
normal  distribution,  S  using  a  normal  shape  function,  F  as 

1  _m  i 

S(i,j)  =  ~Y  Y  ^(i  + 1,  j  +  m) ■  F(l,m) 

C  j=-m  i=~l 

where  F(Z,m)=  Q~8'r  (3) 

m  l 

c=  Y 

j--m  i=-l 

r(Z,  m )2  =  [x(l,  m)  -  x(0, 0)]2  +  [y(l,  m)  -  y(0, 0)]2 


The  newly  mapped  scalar  field,  S  (convoluted  flow  field)  is  computed  by  superposition  of  a  normal  shape 
function,  F,  over  a  sub-grid  consisting  of  2Z+1  by  2m+l  grid  points.  Typically,  a  value  of  7  is  used  for  Z 
and  m,  and  the  grid  point  (0,  0)  denotes  the  peak  location  of  the  shape  function,  F.  The  factor  c  is  a  scale 
factor  to  normalize  the  convoluted  flow  field.  For  example,  the  convoluted  flow  field  will  result  in  a  value 
of  unity  if  X2  =  1  everywhere  in  the  sub-grid  enclosed  by  2/+1  by  2m+l  grid  points.  The  factor  g,  which  is 
a  shape  factor  to  F  for  a  bell-shaped  distribution,  depends  on  the  non-dimensional  grid  resolution,  Ax  or 
Ay.  In  this  report,  g  =  (Ax/ 0.003)2  is  used.  Note  that  convolution  can  be  taken  with  any  of  the  first  or 
second  order  tensors  instead  of  X2  in  the  above  equation. 

The  vortex  center  of  an  individual  velocity  field  can  be  identified  using  various  vortex 
identification  methods.  A  weighted  mean  of  velocity  gradient  tensors  is  computed  for  each  individual 
image.  A  weighted  mean  of  vorticity  is  called  the  vorticity  — mtroid”,  and  is  defined  as 


y=2AlI±i 

IX,- 


(4) 


A  summation  is  performed  over  the  sub-grid,  which  is  a  small  neighboring  region  of  a  peak  vorticity. 
Typically,  we  define  this  sub-grid  as  the  regions  including  the  top  10%  of  the  peak  vorticity.  A  vortex 
center  using  the  second-order  velocity  gradient  tensors  or  a  convolution  of  velocity  gradient  tensors  can 
be  computed  in  a  similar  manner. 
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After  the  vortex  center  is  computed  for  each  of  the  PIV  images,  we  select  80-95  images  for 
averaging  in  order  to  remove  outliers.  This  corresponds  to  a  Z  score  of  1.28-1.96,  which  was 
recommended  by  HART  II  partners  based  on  their  collective  experience.  The  vortices  with  their  centers 
being  outside  of  a  user-defined  Z  score  are  considered  geometric  outliers.  Alternatively,  the  vortices  with 
their  peak  value  of  vorticity,  Q  or  X2  outside  of  the  specified  threshold  are  considered  strength  outliers. 

Considering  a  normally  distributed  sample,  the  standard  deviation  of  vortex  center  is  defined  in  the 
x-  and  y-directions  as 


'lu-y)2 

AM 


(5) 


The  Xj  and  yt  are  the  vortex  center  in  the  x-  and  y-directions  for  an  individual  PIV  image,  N  is  the  sample 
size  and  the  bar  indicates  the  mean  value.  A  distribution  of  the  sample  data  can  be  eccentric  in  either  the 
x-  or  y-direction.  Therefore,  a  circular  standard  deviation  (i .e.,crr  =  JcP  +  o\  )  may  not  be  sufficient  to 

standardize  the  two  dimensional,  normal  distribution.  So,  we  define  here  a  circular  Z  score  for  an 
isotropic  scaling. 

y  =  z  =  4xT+Y 2  (6) 

ax  ay 

For  conditional  averaging,  each  selected  PIV  image  is  collocated  with  its  vortex  center  as  shown  in  Fig.  5, 
and  a  new  vortex  center  is  computed  only  with  the  selected  80-95  individual  PIV  images.  The  conditional 
averaging  is  a  mandatory  process  for  the  analysis  of  vortex  properties  in  the  HART  II  PIV  database  [11]. 


ROTATION  OF  AVERAGED  VELOCITY  FIELD 

Tung  et  al.  [15]  presented  calculations  of  vortex  core  velocities  using  the  vortex  trajectory  method. 
Burley  et  al.  [16]  presented  two  different  methods  to  determine  the  orientation  of  the  vortex  core  with 
respect  to  the  PIV  measurement  plane.  Specifically,  one  method  uses  the  axial  flow  within  the  vortex  core 
and  the  other  uses  the  tangential  (or  spin)  velocity  within  the  vortex  core  to  define  the  orientation.  In  this 
report,  a  simple  geometry -based  vortex  trajectory  method  is  used  to  determine  the  convective  velocity  and 
the  orientation  of  vortex  core  axis. 

In  the  vortex  trajectory  method,  an  averaged  velocity  field  is  rotated  using  a  transformation  matrix. 
The  apparent  vortex  trajectory  is  constructed  using  the  vortex  centers  derived  from  the  averaged  measured 
velocity  fields,  based  on  the  assumption  that  a  vortex  trajectory  is  normal  to  the  vortex  core  axis  or  the 
PIV  measurement  plane.  This  vortex  trajectory  method  was  employed  by  AFDD  using  the  kinematics  of 
Euler  rotation  by  means  of  a  vector  dyadic  [15].  The  method  computes  Euler  rotation  angles  of  the  vortex 
axis,  taking  into  account  the  effect  of  deformation  of  the  vortex  trajectory.  Later,  an  alternative  approach 
was  used  by  DLR  using  the  direction  cosines  computed  from  the  estimated  vortex  trajectory  polynomials 
that  were  constructed  from  neighboring  PIV  measurement  plane  origins  [17]. 

AFDD  Vortex  Trajectory  Method 

For  any  three  adjacent  vortex  positions  Ph  P2  and  P3 ,  one  can  fit  a  circular  arc  through  these  points. 
These  vortex  positions  are  the  vortex  centers  derived  from  the  averaged  velocity  fields.  In  Fig.  6,  the 
normal  vector  (B3)  is  defined  as  the  unit  vector  parallel  to  the  vector  connecting  the  two  positions  (Pj  and 
P3 ),  crossing  at  the  position  P2.  This  normal  vector  and  the  in-plane  vector  connecting  the  positions  Pj 

and  P2  ( PXP2 )  form  the  vortex  trajectory  plane,  and  the  cross  product  of  these  two  vectors  gives  the 
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vertical  vector  (B2)  at  the  position  P2.  Again,  the  cross  product  of  the  vertical  vector  and  the  normal 
vector  at  the  position  P2  gives  a  bi-normal  vector  (Bi).  Thus,  the  rotated  PIV  plane,  which  is  normal  to  the 
vortex  trajectory,  contains  the  bi-normal  vector  in  the  x  direction  (Bi)  and  the  vertical  vector  in  the  y 
direction  (B2).  The  mathematical  expressions  for  these  three  vectors  are  described  in  terms  of  the 
positions,  Pl9  P2  and  P3  as: 


PP 

b3=^,b2  =  b3 

^1^3 


PP 
1  \l  2 

PP 

r\r2 


®1  "  ®2 


B, 


(7) 


Note  that  determination  of  the  normal  vector  (B3)  may  introduce  a  large  error  if  the  position  P2  is  not  the 
mid-point  of  a  circular  arc  connecting  the  three  vortex  positions  (Plt  P2,  P3). 

Rotation  angles  can  be  found  from  kinematics  of  Euler  rotation.  A  transformation  matrix  of  an 
unrotated  PIV  image  to  a  rotated  one  is  determined  from  Euler  rotation  sequence.  The  rotation  angles  are 
computed  using  a  transformation  matrix  in  the  following  dyadic  relation: 


B  =T 


BA 


BA 


■aj=t^ 


A . 


Uj  =12,3  where  =B 


'A  j=BiAj 


(8) 


where  A  and  B  are  the  unit  vectors  in  the  unrotated  and  rotated  planes,  respectively.  The  rotation  tensor 
is  a  transformation  matrix  (dyadic  tensor)  containing  all  three  rotation  angles  in  nature  with  a  rotation 

sequence  about  the  z-,  y-,  and  then  x-axes.  To  eliminate  unnecessary  z-rotation  in  the  rotated  PIV  image, 
one  could  take  the  inverse  z  transformation  such  that 


_  rjnzi/i  rp 

Axy  —  Axyz  z 


i  BA 


(9) 


The  rotation  angle  about  the  z-axis  is  computed  using  T”’1  by  the  relation,  0  -  tan  1  (Tv  /Tn).  Thus,  the 
transformation  matrix  having  a  rotation  about  the  y-axis  and  then  the  x-axis,  is  given  by 


cy 

0 

-sy 

sx-  sy 

cx 

sx-cy 

cx  •  yy 

-sx 

cx'cy 

(10) 


The  finite  rotation  angles  (horizontal  and  vertical)  are  computed  by 


6,  =tan 


-*-23 

V^33  J 


o,=- 


-sm 


do 


This  method  requires  three  vortex  measurement  positions  to  compute  the  associated  vectors.  When 
actual  measurement  is  not  made  at  three  points  along  the  vortex  trajectory,  there  is  not  enough 
information  to  compute  the  vortex  trajectory.  In  such  cases,  the  vortex  positions  could  be  estimated  using 
the  available  PIV  measurement  plane  positions,  ignoring  deformation  of  the  vortex  trajectories. 

DLR  Vortex  Trajectory  Method 

This  method  estimates  the  vortex  trajectory  by  constructing  polynomials  based  on  PIV 
measurement  plane  positions,  based  on  the  assumed  vortex  centers  as  the  origin  of  the  PIV  measurement 
plane  [17].  The  direction  cosines  are  computed  from  the  estimated  vortex  trajectory  polynomials  in  the 
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wind  tunnel  coordinate  system,  and  the  corresponding  direction  cosine  angles  are  inversely  computed  by 
the  following  relations: 


Yx.hub=  tan" 


[  Ay 


hub 


v-A xhubj 


.  and  ry,huh  =  tan1 


-Ax, 


hub 


v  y 


(12) 


These  angles  are  transformed  back  into  the  PIV  measurement  plane  coordinate  system,  where  the  sense  of 
rotation  is  defined  as  follows: 


n  \-rx,h«b’  V  advancing  side 

iV  =  <!  , 

Jxj*b>  V  retreating  side 

Note  that  these  angles  are  direction  cosine  angles  and  can  be  converted  into  Euler  angles  using  a 
coordinate  transformation  matrix. 

The  AFDD  and  DLR  vortex  trajectory  methods  provide  a  good  estimate  for  rotation  angles  of  a 
PIV  image,  and  quantitative  comparison  will  be  given  in  a  later  section.  There  are  the  cases  having  the 
unbalanced  (non-concentric)  contours  of  the  cross  flow  field  surrounding  the  vortex  center,  which  may 
require  an  additional  adjustment  in  the  rotation  of  the  velocity  field.  In  such  cases,  we  incrementally 
rotate  the  unbalanced  cross  flow  field  until  it  becomes  concentric.  When  the  cross  flow  field  is  concentric 
around  the  vortex  center,  the  vorticity  field  also  becomes  concentric  under  the  ideal  situation  [11].  In 
reality,  interactions  with  shear  layers,  turbulence  or  other  unsteady  behaviors  will  make  the  contours  to  be 
non-concentric. 


~Yy,hub  +30°>  V  advancing  side 
Yyhub  -30°>  V  retreating  side 


COMPUTATION  OF  GRADIENTS 

All  flow  field  tensors  such  as  vorticity,  Q  and  are  based  on  the  gradients  of  the  velocity  field,  i.e. 
0V/0x,  9U/dy,  dY/dy  and  dU/dx.  Various  numerical  schemes  are  available  to  compute  these  gradients. 
The  most  commonly  accepted  scheme  is  the  2nd-order-accurate  central  difference  scheme  in  one 
dimension.  More  accurate  schemes  are  available,  such  as  the  3rd-order-accurate  Richardson  extrapolation. 
These  schemes  are  given  as  follows: 


central  difference:  ^ 

dx 

dV 

Richardson  extrapolation:  — 
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(14) 


An  alternative  numerical  scheme  is  to  include  two-dimensional  dependency  of  the  velocity  fields  in 
the  PIV  measurement  planes.  Reference  [18]  provides  a  derivation  based  on  circulation  around  a  point  of 
interest.  The  vorticity  is  expressed  using  Stokes4  theorem  as  follows: 
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z,  1J 
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where  Ay  is  the  area  enclosed  by  the  grid  points,  ( ),  ( i+l,j-l ),  ( i+l,j+l ),  and  (i-JJ+J)  around  the 
point  of  interest.  This  integration  takes  into  account  the  velocity  data  from  eight  surrounding  grid  points. 
For  equally  spaced  grid  points  around  the  point  (/',/).  the  gradients  o'V/fx  and  <5U/<3y  are: 


dV_ 

dx 

dU 

dy 


y  +2V  +V  -V  -2V  -V 


8  Ax 
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Note  that  U  and  V  are  the  velocities  in  the  x-  and  y-directions,  respectively.  This  approach  is 
equivalent  to  applying  the  central  difference  scheme  to  a  smoothed  3x3  kernel  velocity  field,  and  the 
uncertainty  in  the  vorticity  estimate  with  this  approach  is  therefore  expected  to  decrease.  Other  gradients 
such  as  shear  strain  can  be  computed  in  a  similar  fashion.  It  is  noted  that  the  robustness  of  this  approach 
over  the  central  difference  or  Richardson6  s  extrapolation  is  described  in  Ref.  [18]. 


VATISTAS  MODEL 

Tip  vortices  often  interact  with  other  flow  structures  including  shear  layers,  turbulence  and  vortices 
shed  from  other  blades.  A  characterization  of  a  vortex  core  model  can  be  made  with  a  curve-fit  to  a  swirl 
(tangential)  velocity  profile  of  a  vortex  [11,  19].  A  simple  model  is  the  Rankine  vortex,  which  models  the 
core  as  a  solid  body.  The  Oseen-Lamb  model  is  expressed  in  terms  of  the  exponential  decay  function  as  a 
simplified  solution  of  the  Navier-Stokes  equation  [20]: 

r  1  -  e~ap 2  r 

Ve  =  — - - where  p  =  —  (17) 

2  nrc  p  rc 


Tv  is  the  vortex  circulation,  and  rc  is  the  vortex  core  radius.  The  coefficient  a  depends  on  the  kinematic 
viscosity  and  vortex  age,  which  can  be  set  as  1.25643  [19].  A  more  general  mathematical  expression  for  a 
swirl  velocity  profile  is  given  by  Vatistas  et  al.  [21]  as 


2;rre  (1  +  p2n)1,n 


(18) 


At  the  vortex  core  radius  (p=  1),  the  maximum  swirl  velocity  should  satisfy  the  relation, 

r. 


v  =  _ul_  .  2 

v  6>,max  ^  ^ 

2  nr 


-1  In 


(19) 


Therefore,  the  total  circulation  Too  of  a  potential  vortex  is  determined  such  that 


:2^rc’Ve,  max  =V-2 


-1  In 


(20) 


Note  that  the  Scully  core  model,  which  is  widely  used  in  the  rotorcraft  field,  is  identical  to  the  Vatistas 
model  with  n  =  1 .  The  Rankine  vortex  model  is  obtained  as  n  — >  oo,  and  the  Lamb-Oseen  model  is  closely 
approximated  when  n  =  2  [19]. 
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RESULTS  AND  DISCUSSION 


The  PIV  analysis  results  collected  from  HART  II  partners  were  presented  by  the  DLR  [22].  Sample 
comparisons  of  the  derived  vortex  core  properties  for  the  baseline  case  are  presented  in  Fig.  7.  Vortex 
positions,  rotation  angles,  peak  vorticity,  core  radius  and  peak  swirl  velocity  are  shown.  Comparisons  are 
made  in  the  longitudinal  cutting  plane  of  y/R  =  0.7  on  the  advancing  side.  Small-to-large  differences 
among  the  partners6  results  were  found,  and  assessing  the  physical  significance  of  these  differences  was 
difficult.  In  fact,  the  findings  were  not  sufficient  to  resolve  the  differences.  This  report  will  discuss  and 
summarize  the  parameters  and  techniques  which  could  contribute  to  these  differences. 

Discussion  about  the  simple  averaging  technique  [11,  13]  is  not  included  in  this  report  since  HART 
II  partners  have  agreed  that  conditional  averaging  is  required  for  HART  II  PIV  analysis.  The  HART  II 
partners  have  also  agreed  to  use  the  PIV  measured  data  in  24  pixels  x  24  pixels  (interrogation  window 
size  with  50%  overlap)  due  to  more  availability  of  the  data  in  the  HART  II  database.  In  this  report,  only 
the  data  from  the  DLR  measurement  window  will  be  presented.  The  cameras  for  a  close-up  view  (DLR 
window)  had  a  resolution  of  1280  x  1024  pixels,  the  exposure  time  of  an  individual  PIV  image  was  about 
120  ps,  and  the  time  delay  between  the  two  PIV  images  was  mainly  17  ps. 

Conditional  Averaging 

At  each  PIV  measurement  position,  one  hundred  individual  images  are  available.  A  young  vortex, 
which  sheds  from  the  preceding  blade  at  an  azimuth  of  135.6°  (refer  to  Fig.  2a),  quickly  reaches  position 
17  (Dpt  1003)  at  a  wake  age  of  25.3°,  and  is  accompanied  by  a  strong  shear  layer  generated  along  the 
blade  trailing  edge.  For  a  lightly  loaded  rotor,  the  vorticity  strength  of  a  shear  layer  has  a  similar 
magnitude  to  the  tip  vortex  as  in  this  case.  For  a  highly  loaded  rotor  as  seen  on  the  retreating  side  or  in 
hover,  the  vorticity  strength  of  a  shear  layer  becomes  less  significant  compared  to  the  tip  vortex. 

After  passing  the  position  17,  the  young  vortex  continuously  travels  downstream.  When  it  reaches 
the  first  quadrant  of  the  rotor  disk,  the  vortex  is  already  grown  to  a  fully-developed  vortex  and  encounters 
significant  blade -vortex  interactions  (BVI)  with  following  blades,  which  is  a  typical  phenomenon  in 
descending  flight.  A  good  example  of  a  fully  grown  vortex  is  the  one  at  position  22  (Dpt  998),  which  is  at 
a  wake  age  of  425.3°.  Figure  8  shows  typical  (unrotated)  vorticity  maps  for  positions  17  and  22.  The 
vorticity  was  non-dimensionalzed  with  the  rotor  operating  speed.  A  fully  grown  vortex  at  the  position  22 
displays  its  strength  more  diffused  than  a  young  vortex  at  the  position  17. 

The  vortex  centers  of  one  hundred  individual  images  are  plotted  in  Fig.  9  for  position  17  (Dpt  1003) 
of  the  baseline  case.  AFDD  used  a  centroid  for  the  center  identification  while  DLR  used  convolution  of 
X2.  The  centroid  method  displays  individual  vortex  centers  largely  collocated  with  some  scatter,  while  the 
convolution  of  X2  method  results  in  much  less  scatter. 

Figure  10  compares  conditionally  averaged  vortex  centers  for  position  17  (Dpt  1003).  The  vortex 
centers  computed  using  a  centroid  (AFDD),  Q  criterion  (Onera),  and  X2  and  convolution  of  X2  (DLR)  are 
almost  collocated,  with  difference  of  0.02-0.04%  R  (0.36-0.61%  c).  A  tabulated  comparison  of  the  vortex 
properties  for  position  17  is  given  in  Table  la,  which  includes  the  maximum  vorticity,  vortex  centers  and 
their  standard  deviation.  Note  that  1%  of  the  blade  radius  is  equivalent  to  16.5%  of  the  blade  chord 
length. 

Differences  in  vortex  properties  among  various  methods  are  given  in  Table  lb.  The  reference  value 
was  taken  from  the  centroid  result.  The  differences  of  maximum  vorticities  fall  within  7%  of  the  rotor 
speed.  The  vortex  centers  that  were  derived  using  a  centroid,  Q  criterion  and  X2  show  good  agreement,  but 
the  vortex  center  from  convolution  of  X2  is  different  from  the  other  methods  by  about  0.3-0. 5%  of  the 
chord. 

The  vorticity  maps,  computed  using  a  centroid,  Q  criterion,  X2  and  convolution  of  X2,  are  shown  in 
Fig.  11,  exhibiting  a  strong  shear  layer.  This  shear  layer,  resulting  from  airfoil  thickness,  includes  the 
effects  of  the  boundary  layer  (mainly  in  the  cross-flow  component)  as  well  as  the  radial  circulation 
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distribution  (vorticity).  The  vorticity  and  shear  layer  contours  of  the  four  methods  appear  almost  identical 
to  each  other,  except  for  the  value  of  the  maximum  vorticity.  Figure  12  shows  the  same  comparison  but 
with  the  first-  and  second-order  velocity  gradient  tensors,  and  demonstrates  that  the  level  of  difficulty  in 
visual  identification  of  vortex  center  significantly  reduces  with  the  second-order  tensors.  Note  that  the 
coordinate  system  of  the  contour  using  Q  criterion  was  rearranged  at  the  vortex  center.  Although  all  the 
second  order  tensors  (Q  criterion,  Z2  and  convolution  of  Z2)  display  a  cleaner  image  of  the  vortex  core 
than  the  first-order  tensor,  the  accuracy  in  the  vortex  center  location  does  not  depend  on  the  order  of  the 
tensor  —  at  least  within  this  example  —  as  seen  in  Table  1 . 

The  conditionally  averaged  vortex  centers  at  position  22  are  compared  in  Table  2.  The  centroid 
results  are  used  as  a  reference.  The  vortex  center  using  Q  criterion  differs  largely  by  2. 0-3. 3%  of  the 
chord.  Considering  vortex  core  radius  to  be  4-15%  of  the  chord  for  a  typical  rotor  tip  vortex,  this 
difference  is  not  trivial. 

Number  of  Individual  PIV  Images 

When  inspecting  each  individual  PIV  images,  some  images  were  found  to  be  too  diffused  or  blurry. 
In  most  cases,  these  images  were  identified  outliers.  A  typical  characteristic  of  these  outlier  images  is  that 
the  vortex  core  image  is  diffused  and  the  vortex  center  is  away  from  the  mean  vortex  center.  Figure  9 
shows  the  vortex  centers  of  all  100  individual  images  using  the  centroid  and  the  convolution  of  X2.  Outlier 
images  were  eliminated  in  conditional  averaging  (geometric  outlier  removal).  Selecting  good  images  from 
the  population  of  100  images  is  controlled  by  a  Z  score  of  the  normally  distributed  vortex  centers. 
Choosing  the  images  with  a  95%  confidence  level  corresponds  to  a  Z  score  of  1.96,  and  HART  II  partners 
agreed  to  use  about  80-95  images,  corresponding  to  Z  =  1.28-1.96. 

Table  3  shows  a  relationship  between  a  Z  score  and  the  number  of  the  PIV  images  selected.  As 
expected,  the  maximum  vorticity  rapidly  decreases  when  more  images  are  included  for  averaging.  The 
corresponding  vorticity  maps  are  shown  for  position  17  in  Fig.  13.  The  vortex  contours  are  almost 
unchanged  with  a  variation  of  Z  scores,  but  the  peak  vorticity  values  (hot  spot)  are  moderately  changed. 
With  39  images  (Z=0.5),  the  peak  vorticity  increases  by  16%  from  the  vorticity  with  89  images. 

In  addition  to  this  geometric-position-based  elimination  of  outliers,  one  can  alternatively  consider 
dealing  with  vortex  strength  (strength  outlier  removal)  by  analyzing  the  peak  vorticity,  peak  Q,  peak  Z2  or 
the  peak  value  of  the  convolution  with  any  of  those  inside  a  vortex  core.  For  example,  one  can  eliminate 
the  images  having  unusually  lower  vortex  strength  by  defining  a  threshold  (or  Z  score)  below  the  mean. 
The  reasoning  is  that  blade-vortex-interaction  noise  is  dominated  by  the  strongest  interactions  and  a 
proper  averaging  should  thus  retain  the  strongest  vortices  for  analysis  of  their  properties.  The  number  of 
the  images  eliminated  by  this  method  is  similar  to  the  number  with  the  geometric  outlier  removal,  and 
thus  the  number  of  remaining  images  with  the  strength  outlier  removal  is  typically  in  the  range  of  80-95. 
However,  in  cases  where  the  vortex  of  interest  is  close  to  another  one  having  the  opposite  sense  of 
rotation  (dual  vortices),  the  number  of  images  may  drop  down  significantly,  and  an  automatic  procedure 
for  robust  identification  of  the  images  requires  more  complicated  sub-processes. 

Sub-grid  in  Vortex  Center  Identification 

A  vorticity  field  typically  displays  a  hot  spot  where  a  strong  vortex  lies,  and  so  it  is  not  difficult  to 
visually  identify  a  vortex  center.  An  automatic  identification  of  a  vortex  center  using  a  computer  program, 
though  more  consistent  and  objective,  is  not  as  straightforward  as  a  visual  identification.  It  is  more 
difficult  to  identify  a  vortex  center  containing  some  noise  and  a  shear  layer.  Figure  14  shows  a 
perspective  view  of  the  vorticity  map  of  image  #3  at  position  17.  A  sharp  peak  of  a  vortex  is  accompanied 
by  a  series  of  moderate  peaks  from  a  shear  layer. 

The  PIV  measurement  data  are  provided  on  a  grid,  and  the  scalar  values  (first  and  second  order 
tensors)  are  computed  on  the  grid  to  identify  a  vortex  center.  The  physical  vortex  center,  however,  is  not 
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confined  to  any  resolution  represented  by  this  grid,  and  can  be  anywhere  between  the  grid  points.  For  this 
reason,  re-meshing  of  the  grid  points  by  interpolation  is  required. 

A  vortex  center  in  the  vorticity  map  is  determined  by  first  computing  the  vortex  strength  (yz,  X2,  Q 
or  convolution  of  X2)  on  the  grid  and  then  searching  for  the  maximum  vortex  strength  to  determine  the 
vortex  center.  A  sub-grid  is  a  small  neighboring  region  surrounding  peak  vorticities,  and  the  neighboring 
region  is  defined  by  a  specified  threshold  in  terms  of  the  maximum  vortex  strength.  Figure  14  shows  sub¬ 
grids  with  the  top  10%  and  40%  of  the  vorticity  peak,  and  the  latter  case  (top  40%)  displays  two  sub-grids 
which  would  result  in  a  large  shift  of  the  vortex  center. 

Figure  15  shows  the  vortex  centers  of  image  #3  at  position  17  with  sub-grids  defined  using  the  top 
10%  to  50%  of  the  peak  vorticity.  Note  that  the  grid  resolution  in  the  PIV  data  was  1.25%  of  the  chord  (= 
0.0756%  R)  in  both  the  x-  and  y-directions.  The  sub-grid  with  the  top  50%  results  in  a  vortex  center 
shifted  15%  chord  away  from  the  vortex  center  computed  with  the  top  10%.  As  observed,  the  sub-grid 
size  is  sensitive  in  determining  a  vortex  center.  The  question  remains  about  what  sub-grid  size  is  adequate 
for  determining  the  vortex  center.  We  recommend  using  the  top  10%  of  the  maximum  strength  of  a 
vortex  (yz,  X2,  Q  or  convolution  of  X2)  for  a  sub-grid  definition.  Note  that  this  sub-grid  is  used  only  for 
identifying  the  vortex  centers  in  each  individual  image. 

Figure  16  shows  the  perspective  and  top  views  of  X2  and  convolution  of  X2  of  image  #3  at  position 
17  (Dpt  1003).  As  compared  to  the  vorticity  using  a  centroid  method  in  Fig.  14,  the  perspective  view  of  X2 
(Fig.  16a)  seems  slightly  cleaner.  But,  the  views  of  the  convolution  of  X2  display  less  noisy  even  when 
compared  with  the  same  2nd-order  tensor,  X2. 


Conditional  Averaging  Using  a  Subset  of  the  PIV  Database 

The  PIV  data  analysis  is  an  integral  of  all  the  sub-level  processes  in  conditional  averaging  with  the 
various  numerical  differentiation  schemes  for  computing  velocity  gradients,  which  affect  the  first-  and 
second-order  tensor  data.  Establishing  thresholds  for  removing  outliers  is  another  sub-level  process 
followed  by  applying  various  vortex  center  identification  methods.  The  differences  between  HART  II 
partners6  results  (Fig.  7)  originated  from  each  of  the  sub-level  processes.  Due  to  the  large  size  of  the  PIV 
database,  data  analysis  became  very  complex  and  quantifying  these  differences  at  each  sub-process  level 
became  challenging.  To  help  this  situation,  we  attempted  a  conditional  averaging  with  a  subset  of  the 
data,  so  that  a  few  PIV  images  were  used  instead  of  all  100  images. 

For  this  exercise,  two  PIV  positions  were  selected  in  the  longitudinal  plane  of  y/R  =  0.7  on  the 
advancing  side:  position  17  (Dpt  1003)  and  position  22  (Dpt  998).  Two  individual  images  were  used  for 
position  17  and  three  images  for  position  22.  Conditionally  averaged  results  with  these  images  are 
compared  for  the  two  positions  17  and  22. 

a)  Case  Is  Pos  17  (Dpt  1003) 

Figure  17  shows  the  first  16  individual  PIV  images  at  position  17  of  the  baseline.  Among  these  16 
images,  images  #1  and  #3  are  selected  for  this  exercise.  Figure  18  shows  a  comparison  of  vortex 
properties  for  image  #1  as  calculated  by  the  HART  II  partners.  The  figure  shows  a  strong  vortex  (hot 
spot),  accompanied  by  a  moderately  strong  shear  layer.  The  maximum  vorticity  is  much  higher  than  the 
conditionally  averaged  value  using  all  the  images  (see  Table  1).  The  label,  in  the  contour  stands  for 
AFDD,  J36  for  DLR,  and  =06  for  Onera.  In  a  large  window  (left  image),  vortex  centers  are  collocated 
almost  to  each  other,  and  so  this  image  is  enlarged  by  roughly  25  times  (right  image).  Note  that  one  grid 
spacing  in  this  enlarged  image  is  equal  to  1.25%  of  the  chord  (=  0.0756%  R),  and  1%  of  the  radius  (%  R) 
equals  to  16.5%  of  the  chord.  The  maximum  vorticity  difference  (Amax)  between  different  methods  is 
small,  and  the  difference  in  the  circular  vortex  center  (Arv  =^Ay2  +  A y]  )  is  0.06%  R  (=  0.99%  c),  as 
shown  in  Table  4a. 
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A  similar  comparison  was  made  for  image  #3  in  Fig.  19  and  Table  4b.  This  image  appears  slightly 
fuzzy.  Computed  vortex  properties,  such  as  vorticity  as  well  as  vortex  center,  are  almost  identical  to  each 
other  when  using  the  centroid  and  Q  criterion  methods.  The  vorticity  by  the  convolution  of  X2  is  also 
similar  to  the  values  by  other  methods,  but  its  vortex  center  shifts  to  the  x-direction  by  about  0.08%  R  (= 
1.36%  c). 

Finally,  both  images  were  used  for  conditional  averaging.  The  vortex  properties  with  images  #1  and 
#3  are  compared  in  Fig.  20  and  Table  4c-4d.  The  vorticity  by  the  convolution  of  X2  is  10%  higher  than  the 
value  using  a  centroid  or  Q  criterion.  The  vortex  center  by  the  convolution  of  X2  differs  from  the  centroid 
method  result  by  0.03%  R  (=  0.50%  c). 

b)  Case  2:  Pos  22  (Dpt  998) 

Figure  21  shows  first  16  individual  images  at  position  22  of  the  baseline.  Images  #2,  #3  and  #8  are 
used  for  conditional  averaging.  The  images  #2  and  #3  show  a  clear  hot  spot  and  the  image  #8  seems 
slightly  fuzzy.  The  vortex  at  this  position  is  fully  grown,  its  vortex  strength  diffused,  and  the  flow  field 
seems  relatively  noisy. 

Vortex  properties  are  derived  using  conditional  averaging  by  the  centroid  and  convolution  of  X2 
methods  shown  in  Fig.  22  and  Table  5.  The  peak  vorticity  at  the  position  22  is  substantially  less  than  the 
peak  at  the  position  17.  The  circular  vortex  center  is  different  by  0.054%  R  (=  0.89%  c)  between  the  two 
methods.  Note  that  a  fuzziness  of  the  vortex  image  in  this  figure  is  seen  typically  from  older  vortices. 

Rotation  Angles  for  an  Averaged  Image 

The  PIV  measurement  plane  is  not  generally  normal  to  the  vortex  trajectory  due  to  pre-setting  the 
orientation  of  the  measurement  plane  prior  to  the  test.  Since  an  accurate  analysis  of  the  vortex  core  can  be 
made  only  in  the  vortex  axis,  the  rotation  angles  from  the  measurement  plane  to  the  vortex  axis  should  be 
accurately  identified.  In  this  study,  these  rotation  angles  were  defined  in  terms  of  Euler  angles  with  a 
sequence  of  the  y-x  rotation,  or  direction  cosine  angles.  Ideally,  the  resulting  vortex  image  will  be  normal 
to  the  vortex  trajectory  after  rotation. 

As  a  general  hypothesis,  every  vortex  in  an  individual  PIV  image  is  assumed  to  have  the  same 
inclination  of  the  vortex  axis  with  respect  to  the  measurement  plane,  though  the  vortex  axis  orientation 
may  vary  in  each  individual  image.  During  the  PIV  analysis  exercise  among  the  HART  II  partners,  the 
discussion  evolved  whether  individual  images  should  be  rotated  prior  to  conditional  averaging  or  not,  and 
the  conclusion  was  that  this  would  introduce  additional  complication  and  uncertainties  without  having 
significant  benefits.  After  the  rotation  of  the  velocity  field,  the  meshes  in  the  rotated  image  are  no  longer 
orthogonal  and  thus  data  interpolation  on  the  rotated  grid  of  the  individual  images  becomes  more 
complicated.  In  addition,  the  analysis  of  rotation  angles  for  each  individual  is  repeated  over  the  number  of 
individual  images.  Reference  [11]  showed,  using  hover  data,  that  averaging  of  individual  rotation  angles 
was  in  agreement  within  one  degree  tolerance  despite  a  standard  deviation  of  about  four  degrees,  and  this 
difference  is  considered  negligible.  Therefore,  the  conditionally  averaged  velocity  field  will  be  used  in  the 
present  work. 

Figure  23  shows  the  contours  of  vorticity  (y/Q)  and  cross  flow  (Vzv/QR)  before  and  after  rotation. 
Note  that  the  cross  flow  velocity  was  filtered  out  by  the  velocity  at  the  vortex  center.  The  conditionally 
averaged  image  was  generated  using  a  centroid  method  with  images  #2,  #3  and  #8  at  position  22.  The 
rotation  angles  were  found  to  6.3°  and  -34.3°,  respectively,  in  the  x-  and  y-axes  of  the  PIV  measurement 
plane.  The  peak  vorticity  of  the  unrotated  image  was  24.40,  and  that  of  the  rotated  image  was  24.41. 
Small  differences  are  observed  in  both  peak  values  and  the  contour  patterns,  which  forces  to  preclude 
using  the  vorticity  map  forjudging  the  quality  of  rotation  angles.  On  the  contrary,  the  cross  flow  velocity 
field  shows  a  relatively  good  quality  of  concentric  contour  after  rotation.  Therefore,  the  cross-flow 
velocity  field  is  considered  a  primary  indicator  for  an  identification  of  rotation  angles  in  this  report. 
Though  the  resulting  image  is  a  projection  of  the  measured  data  at  a  view  angle  parallel  to  the  vortex  axis, 
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it  is  assumed  that  the  flow  field  may  not  be  significantly  different  within  the  volume  covered  by  the  z- 
coordinate  of  the  image  after  rotation.  The  projection  can  be  also  assumed  as  being  measured  in  the  plane 
normal  to  the  vortex  axis,  and  thus  an  analysis  of  vortex  properties  may  be  allowed  within  this  projected 
data  field. 

Figure  24  compares  two  rotation  angles  (horizontal  and  vertical),  at  position  22,  derived  using  the 
centroid  (AFDD)  and  the  convolution  of  X2  (DLR).  Larger  differences  in  the  rotation  angles  between  the 
two  methods  are  found  from  image  #8.  This  result  is  expected  since  image  #8  is  fuzzy  compared  with  the 
other  two  images  and  so  more  difficult  to  accurately  identify  the  vortex  properties.  The  rotation  angles 
after  being  conditionally  averaged  (ca)  are  shown  in  the  figure.  Interestingly  the  difference  of  the  rotation 
angles  from  the  averaged  image  is  less  than  that  from  the  single  image  #8,  and  a  similar  example  is  found 
in  Ref.  [11].  In  Table  6,  the  difference  in  the  rotation  angles  between  the  two  methods  is  about  3°  or  less 
with  a  conditional  average,  but  the  difference  increases  to  6°  to  12°  when  the  rotation  angles  are  compared 
between  individual  images.  Derived  rotation  angles  using  all  100  images  (see  Fig.  7),  are  given  in  Table  7 
for  positions  17  and  22.  Note  that  after  rotation  into  the  vortex  axis  system,  the  resulting  grid  is  no  longer 
orthogonal  and  skewed  instead.  Care  must  be  taken  when  derivatives  such  as  dV/dx  are  computed,  since 
the  data  must  be  properly  interpolated. 

Swirl  Velocity 

An  accurately  rotated  vorticity  or  velocity  field  in  the  PIV  image  displays  a  concentric  contour 
around  the  vortex  center,  provided  that  the  vortex  is  rotationally  symmetric.  This  concentric  contour  can 
be  distorted  due  to  interactions  with  other  vortices,  shear  layers,  or  turbulence.  From  the  rotated,  averaged 
velocity  field,  a  swirl  velocity  can  be  reconstructed  in  two  steps.  First,  tangential  velocities  are  collected 
from  the  cutting  axes  —  the  straight  lines  in  the  radial  direction  from  a  vortex  core.  The  direction  of  the 
cutting  axes  can  be  to  the  right,  up,  left,  down  and  in  between,  and  this  process  generates  multiple 
tangential  velocity  profiles.  Secondly,  by  taking  an  algebraic  averaging  of  these  velocity  profiles,  the 
desired  tangential  velocity  profile  is  obtained. 

As  an  alternative  approach,  all  individual  images  can  be  separately  analyzed  instead  of  a  single 
averaged  image.  Collecting  individual  velocity  profiles  provides  a  tangential  velocity  — loud”,  from 
which  an  averaged  tangential  velocity  profile  can  be  computed.  The  outlier  removal  methodology  can  be 
applied  to  reduce  the  bandwidth  of  the  cloud  before  averaging  (DLR  approach). 

Figure  25  shows  the  vorticity  map  and  tangential  velocities  of  the  conditionally  averaged  image  at 
position  17  using  the  centroid  method.  Images  #1  and  #3  were  used  for  averaging.  The  cutting  axis 
azimuths  are  shown  on  the  vorticity  map  as  an  example,  and  the  associated  tangential  (swirl)  velocity 
profiles  are  also  presented.  The  cloud  in  the  tangential  velocity  profiles  represents  all  the  tangential 
velocities  possible  from  the  averaged  image  at  this  position  17,  and  the  tangential  velocities  show  a  large 
scatter  which  is  likely  due  to  a  shear  layer.  Figure  26  shows  a  similar  plot  for  position  22,  where 
averaging  was  performed  with  images  #2,  #3  and  #8. 

To  determine  a  tangential  velocity  profile,  an  appropriate  method  for  choosing  cutting  axis  is 
needed,  and  it  would  take  into  account  the  least  noisy  flow  field  where  the  influence  of  other  neighboring 
vortices  or  shear  layers  is  minimized.  Practically,  a  perfectly  concentric  flow  field  around  the  vortex  core 
is  not  feasible  and  the  judgment  about  the  quality  of  the  vortical  flow  field  is  also  subject  to  a  researchers 
experience.  A  common  approach  for  determining  the  cutting  axis  (horizontal,  vertical,  in  between,  or  a 
combination  of  these  choices)  is  to  choose  the  axes  manually  in  the  region  where  a  cleaner  vorticity  field 
is  observed.  For  example,  the  optimum  cutting  axis  is  selected  as  the  azimuths  of  90°  and  180°  for  the 
positions  17  and  22  as  in  Figs.  25  and  26,  respectively. 

The  corresponding  tabulated  data  of  vortex  core  radii  and  tangential  velocities  are  given  in  Table  8. 
When  the  four  horizontal  and  vertical  axes  (0°,  90°,  180°  and  270°)  are  selected,  the  mean  vortex  core 
radii  are  5.3%  chord  for  the  position  17  and  10.7%  chord  for  the  position  22.  The  use  of  the  optimum 
cutting  axis  at  90°  and  180°  enables  to  substantially  reduce  the  vortex  core  radius  compared  with  the  use 
of  the  four  axes  —  for  a  young  vortex  (position  17)  by  40%  to  3.2%  chord,  but  for  an  older  vortex 


13 


(position  22)  it  unexpectedly  increases  the  core  radius  by  7%.  The  peak  tangential  velocities  with  the 
optimum  cutting  axes  increase  by  less  than  8%  for  both  positions  17  and  22. 

Figure  27  is  an  example  of  the  curve-fit  results  using  a  Vatistas  model  for  positions  17  and  22,  and 
the  selected  cutting  axes  are  90°  and  180°  for  both  positions.  For  position  17,  the  core  radius  was  3.2% 
chord  (=  0.20  %R)  with  the  peak  swirl  velocity  2.9%  of  the  rotor  tip  speed.  The  Vatistas  exponent  of 
n= 0.60  was  best-fit  to  the  tangential  velocity.  For  position  22,  the  core  radius  was  1 1.4%  chord  (=  0.69% 
R)  with  the  peak  swirl  velocity  2.58%  of  the  tip  speed.  The  Vatistas  exponent  was  slightly  increased  to 
n=0.70.  Note  that  the  Vatistas  exponent  may  depend  on  the  choice  of  the  cutting  axes.  This  exponent  n 
can  be  considered  as  a  shape  factor  for  the  tangential  velocity  profile.  Since  the  number  of  data  points 
within  the  peak  velocity  is  relatively  scarce,  Vatistas  exponent  n  may  not  be  highly  accurate. 


CONCLUSIONS 

Despite  the  past  experiences  of  HART  II  partners  who  have  dealt  with  PIV  data  analyses  for  years, 
the  complexity  of  the  PIV  analysis  procedures  became  readily  apparent  in  the  large  HART  II  database. 
Comparisons  of  the  partners4  results  typically  led  to  unresolved  questions  and  ambiguity.  As  a  result,  the 
partners  focused  on  a  much  smaller  sample  of  the  database.  Although  the  use  of  a  subset  of  the  PIV 
database  may  not  be  sufficient  to  establish  comprehensive  guidelines  for  establishing  a  PIV  data 
processing  methodology,  the  effort  in  this  report  certainly  offers  some  insights  on  the  PIV  data  analysis 
techniques,  and  helps  establish  a  strong  foundation  for  an  analysis  methodology.  The  following  are  of  the 
main  findings  from  this  study: 

1)  Vortex  center  identification  method  -  Four  methods  were  used  for  comparison.  The  methods 
were:  centroid  (AFDD,  DLR),  Q  criterion  (ONERA,  AFDD,  DLR),  X2  (DLR,  AFDD)  and  convolution  of 
X2  (DLR).  Although  the  contours  of  all  the  second  order  velocity  gradient  tensors  (Q  criterion,  Z2  and 
convolution  of  Z2)  displayed  a  strong  hot  spot  surrounded  with  a  cleaner  flow  field  compared  to  the  first- 
order  tensor  (vorticity),  the  accuracy  of  the  computed  vortex  center  location  was  not  affected  by  the  order 
of  the  tensor.  The  difference  of  the  vortex  center  locations  among  the  partners  was  1.5  -  1.7%  of  the  chord 
for  position  17,  and  2.0  -  3.3%  for  position  22. 

2)  Number  of  individual  PIV  images  -  Although  the  vortex  contour  was  almost  unchanged  with  a 
variation  of  the  number  of  images,  the  peak  vorticity  value  (hot  spot)  was  moderately  sensitive.  The  peak 
vorticity  for  position  17  increased  by  16%  with  38  images  (Z=0.5)  compared  to  the  value  with  89  images 
(Z=1.5). 

3)  Sub-grid  in  vortex  center  identification  -  The  sub-grid  with  the  top  50%  of  the  peak  vorticity 
showed  the  vortex  center  shift  by  15%  of  the  chord  from  the  reference  value  computed  with  the  top  10%. 
The  sub-grid  size  was  important  in  determining  a  vortex  center,  and  the  top  10%  of  the  maximum 
vorticity  for  a  sub-grid  definition  was  recommended. 

4)  Conditional  averaging  with  a  subset  of  the  PIV  database  -  Since  a  subset  of  the  PIV  images 
significantly  reduced  the  problem  size,  it  was  possible  that  each  PIV  analysis  technique  was  closely 
examined  and  accordingly  the  differences  between  HART  II  partners4  results  were  better  understood. 

5)  Rotation  angles  -  The  difference  of  the  rotation  angles  between  the  averaged  images  computed 
using  the  centroid  and  convolution  of  Z2  was  less  than  3°  for  position  22,  but  the  differences  between 
individual  images  were  interestingly  larger  (6°  to  12°). 

6)  Swirl  velocity  -  The  optimum  cutting  axis  for  extracting  swirl  (tangential)  velocities  was  90° 
and  1 80°  azimuths  for  the  cases  analyzed,  which  reduced  substantially  the  vortex  core  radius  for  a  young 
vortex  (position  17)  by  40%  to  3.2%  chord  when  compared  to  using  four  cutting  axes  (0°,  90°,  180°  and 
270°).  A  Vatistas  model  gave  good  estimates  for  swirl  velocity  profiles. 
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Table  1.  Vortex  properties  of  the  PIV  image  at  position  17 


a)  Conditionally  averaged  properties 


centroid a 

Q  criterion0 

V 

conv  X2  d 

Anax 

Y  max/O 

33.1 

30.8 

32.7 

31.5 

2.3 

xv/R  (%) 

0.41 

0.43 

0.40 

0.44 

0.04 

yv/R  (%) 

0.52 

0.52 

0.52 

0.54 

0.02 

ax/R(%) 

0.60 

* 

0.61 

0.13 

0.48 

av/R  (%) 

0.29 

* 

0.30 

0.11 

0.19 

b)  Difference  in  the  averaged  property,  referenced  from  a  centroid 


centroid a 

Q  criterion0 

conv  X2  d 

Aymax/£2  (%) 

r 

-6.9 

-1.2 

-4.8 

Axv/c  (%) 

r 

0.26 

-0.15 

0.46 

Ayv/c  (%) 

r 

-0.03 

-0.18 

0.28 

a  provided  by  AFDD, 0  provided  by  Onera,  d  provided  by  DLR, r  reference  data, 
no  data  provided 


Table  2.  Vortex  properties  of  the  PIV  image  at  position  22 


a)  Conditionally  averaged  properties 


centroid a 

Q  criterion0 

V 

conv  X2  d 

Anax 

Ymax/^2 

24.8 

25.6 

24.9 

25.2 

0.8 

xv/R  (%) 

-0.08 

0.04 

-0.08 

-0.13 

0.17 

yv/R  (%) 

0.03 

-0.17 

0.04 

0.01 

0.21 

ax/R  (%) 

0.76 

* 

0.76 

0.57 

0.19 

av/R  (%) 

0.46 

* 

0.46 

0.41 

0.05 

b)  Difference  in  the  averaged  property,  referenced  from  a  centroid 


centroid a 

Q  criterion0 

^2a 

conv  X2  d 

Aymax/f2  (%) 

r 

1.6 

0.4 

1.6 

Axv/c  (%) 

r 

1.98 

0.0 

-0.83 

Ayv/c  (%) 

•  1  1  1  A  rTATX 

r 

O  •  1  1  1 

-3.31 

rv  _  d  •  1  1  1 

0.17 

t\t  n  r  n _ 

-0.33 

a  provided  by  AFDD, 0  provided  by  Onera,  d  provided  by  DLR, r  reference  data, 
no  data  provided 
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Table  3.  Z  score,  the  number  of  images  and  the  maximum  vorticity  for  the  PIV  image  at  position  17 
when  the  centroid  method  is  employed 


Z  score 

#  of  images 

Y  max/ 

(y-y  )/  y 

Z=1.5  Z=1.5 

0.5 

38 

38.3 

15.7% 

0.7 

62 

36.5 

10.3% 

1.5 

89 

33.1 

r 

5.0 

99 

31.7 

-4.2% 

r  reference  data 


Table  4.  Vortex  properties  derived  with  the  individual  PIV  images  #1  and  #3  at  position  17 
a)  Image  #1 


centroid a 

Q  criterion0 

conv  Z2  d 

A  max 

Ymax/^2 

51.9 

52.6 

52.6 

0.7 

xv/R  (%) 

0.31 

0.35 

0.29 

0.06 

yv/R  (%) 

0.52 

0.52 

0.53 

0.01 

b)  Image  #3 


centroid a 

Q  criterion 0 

conv  Z2  d 

Amax 

Ymax/^2 

34.1 

34.1 

34.1 

0.0 

xv/R  (%) 

0.35 

0.35 

0.43 

0.08 

yv/R  (%) 

0.75 

0.75 

0.73 

0.02 

c)  Conditional  averaging  with  images  #1  and  #3 


centroid a 

Q  criterion0 

conv  Z2  d 

Amax 

Y  max/^2 

43.3 

43.4 

47.7 

4.4 

xv/R  (%) 

0.33 

* 

0.36 

0.03 

yv/R  (%) 

0.64 

* 

0.64 

0.0 

d)  Difference  in  conditionally  averaged  property  with  images  #1  and  #3,  referenced  from  the 
results  of  a  centroid 


centroid a 

Q  criterion0 

conv  Z2  d 

Ayma x/n  (%) 

r 

0.2 

10.2 

Axv/c  (%) 

r 

* 

0.50 

Ayv/c  (%) 

r 

* 

0.0 

a  provided  by  AFDD, 0  provided  by  Onera,  d  provided  by  DLR, r  reference  data, 
no  data  provided, 


18 


Table  5.  Vortex  properties  derived  with  the  individual  PIV  images  #2,  #3  and  #8  at  position  22 


a)  Conditional  averaging 


centroid a 

conv  X2  d 

^max 

YmaxA2 

24.4 

23.3 

1.1 

xv/R  (%) 

-0.29 

-0.34 

0.05 

yv/R  (%) 

0.35 

0.37 

0.02 

b)  Difference  in  conditionally  averaged  property,  referenced  from  a  centroid 


centroid a 

conv  X2  d 

AYma */n  (%) 

r 

4.5 

Axv/c  (%) 

r 

0.83 

Ayv/c  (%) 

r 

0.33 

a  provided  by  AFDD,  provided  by  DLR, r  reference  data 


Table  6.  Rotation  angles  computed  with  the  individual  PIV  images  #2,  #3  and  #8  at  position  22 
a)  Image-to-image  (i2i)  difference 


centroid a 

conv  X2  d 

A0x,i2i 

7.7° 

12.3° 

A0y,i2i 

6.2° 

5.2° 

b)  Conditional  average  (ca) 


centroid a 

conv  X2  d 

Anax 

0x,ca 

6.3° 

00 

OJ 

o 

2.0° 

0y,ca 

-34.3° 

-31.4° 

2.9° 

a  provided  by  AFDD,  provided  by  DLR 
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Table  7.  Rotation  angles  computed  from  the  conditionally  averaged  image  with  all  100  PIV  images  at 
positions  17  and  22 

a)  Pos  17 


centroid a  Q  criterion 0  conv  X2  d 


b)  Pos  22 


centroid a  Q  criterion 0  conv  X2  d 


a  provided  by  AFDD, 0  provided  by  Onera,  d  provided  by  DLR 


Table  8.  Vortex  core  radii  and  tangential  velocities  at  positions  17  and  22,  averaged  (centroid)  with  a  few 
selected  images 


Pos  17 

Pos  22 

cutting  axis 

rc  /c  (%) 

v0,max/QR  (%) 

rc  /c  (%) 

v0,max/QR  (%) 

0° 

4.8 

2.4 

10.4 

2.1 

o 

o 

2.8 

3.0 

12.4 

2.6 

00 

o 

o 

3.6 

2.9 

10.4 

2.6 

270° 

10.1 

2.6 

9.8 

2.5 

mean 

(0o,90°,180o,270°) 

5.3 

2.7 

10.7 

2.4 

optimum  mean 
(90°,  180°) 

3.2 

2.9 

11.4 

2.6 

20 


b)  Reference  blade  at  20  degrees 

Figure  2.  Schematic  diagram  for  PIV 
measurement  planes 
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Figure  3.  Raw  image  from  two  different  cameras  for  the 
baseline  case,  position  21. 
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Y(%R)  Y(%R) 


Horizontal  Velocity  Map,  U 
Dpt  1003,  Pos  17,  Image  #1 


Vertical  Velocity  Map,  V 
Dpt  1003,  Pos  17,  Image  #1 


Normal  Velocity  Map,  W 
Dpt  1003,  Pos  17,  Image  #1 


Vorticity  Map 

Dpt  1003,  Pos  17,  Image  #1 


Figure  4.  Single  PIV  image  of  3-C  velocity  and  vorticity  map  (image  #1)  for  position  17  (Dpt 
1003)  with  an  age  of  25.3°. 
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100 
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swirl 

velocity 


Figure  5.  Conditional  averaging  procedure 


Figure  6.  Re-orientation  of  the  PIV  measurement  plane  for  the  vortex  trajectory  method. 
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Vortex  Horizontal  Position,  BL,  y/R=0.7 


Vortex  Vertical  Position,  BL,  y/R=0.7 


xwt/R 


xw  t/R 


Horizontal  Angle,  BL,  y/R=0.7 


Vertical  Angle,  BL,  y/R=0.7 


Peak  Vorticity,  BL,  y/R=0.7 


Vortex  Core  Radius,  BL,  y/R=0.7 


Maximum  Swirl  Velocity,  BL,  y/R=0.7 


Figure  7.  Comparison  of  vortex  properties  among  HART  II  partners4  results  in  the 
longitudinal  plane  of  y/R=0.7  on  the  advancing  side  (baseline,  positions  17-23) 
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Y(%R) 


Dpt  1003,  Pos  17,  Centroid  Dpt  998,  Pos  22,  centroid 


Figure  8.  Conditionally  averaged  vorticity  maps  (unrotated)  for  positions  17  (Dpt  1003,  wake  age 
of  25.3°)  and  22  (Dpt  998,  wake  age  of  425.3°). 
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Figure  9.  Vortex  center  distribution  of  all  100  individual  images  for  position  17  (Dpt  1003). 
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Figure  10.  A  comparison  of  the  mean  vortex  center  derived  using  centroid,  Q  criterion  and 
convolution  of  A,?  for  position  17  (Dpt  1003). 
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Figure  11.  Comparison  of  conditionally  averaged  vorticity  (y/Q)  using  various  center 
identification  methods  for  position  17  (Dpt  1003). 
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Figure  12.  Comparison  of  conditionally  averaged  velocity  gradient  tensors  -  centroid,  Q  criterion, 
X2,  and  convolution  of  X2  -  using  various  center  identification  methods  for  position  17  (Dpt 
1003). 
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Figure  13.  Effect  of  Z  scores  or  standard  deviation  of  vortex  center  on  the  conditionally  averaged 
vorticity  for  position  17  (Dpt  1003):  standard  deviations  of  0.5,  0.7,  1.5,  and  5.0. 
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Figure  14.  A  perspective  view  of  the  vorticity  map  and  the  locations  of  its  sub-grids  in  the 
PIV  image  #3  for  position  17  (Dpt  1003)  in  the  baseline  case. 
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17 


a)  Centroid  b)  Q  criterion 

15.  Vortex  centers  with  a  variation  of  the  sub-grid  size  for  the  PIV  image  #3  for  position 
(Dpt  1003). 
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b)  Top  view  of  X2 


d)  Top  view  of  convolution  of  X2 


Figure  16.  Perspective  and  top  views  of  X2  and  convolution  of  X2  in  the  PIV  image  #3 
for  position  17  (Dpt  1003)  in  the  baseline  case. 
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Figure  17.  The  first  16  individual  PIV  images  at  position  17  (Dpt  1003)  in  the  baseline  case. 
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BL,  Pos  17,  image  #1 
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Figure  18.  Individual  Image  #1:  position  17  (Dpt  1003). 
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Figure  19.  Individual  Image  #3:  position  17  (Dpt  1003). 
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Conditional  Averaging  (AFDD) 
BL,  Pos  17,  images  #1,  #3 
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Figure  20. Conditional  averaging  with  images  #1  and  #3:  position  17  (Dpt  1003). 


OjCC 

s  0 


■OjCC 

■001 


•J-  v  it 


%  i  Z 

K&siPni  ■ 

oW  ‘tS.te1  ‘-6b‘*‘  ‘-Ote1  ‘/{t1  ‘dbi 1 M  ‘  '6IA‘  ‘tie 


■0.02 

-0.0* 

•oifcf 


\  V/W.o 

^Sf***-*  •/  »5I  ■'  H 

*  o' 

jefegy  2 1 

?ft:  7^'  •ate'uir 


2  0 
-OjCC 
■OOi 


^  •  r*  .  *7  -  7/, 

-T*  KjJ.j  ^Tvj'  Jr*  002 

;  .  /  ■  *  '  *•  >\  i  o 


' 


^  3s 


•o.oc 

_  ^  Ji 

flLJ  ■■ 

I  jCT'fe  "i'.'tfe'  '0M"6fe'  'His 


Boisys^gg 

K  ££&S3&i 

r  *  3B 
>*c 

4 

^  I 

vri 

-o^1  ‘-dbV  ‘-6fe‘  ‘.‘ ‘  ‘dA*‘  ‘tf  Ai‘  ‘tf  jv*‘  ‘t>h 


-  * ' 


o.oc 

!  0 
•0.02 
-0.0* 
-oSr 


. , «■  , 


6 

‘rfjfeT1  ‘dAl‘  ‘dfe‘  ‘flfe 


?  0 
-o.oc 
-0.0* 


2*'  VCW 


yfa 


1  % 


jrr  >  ^  »£  { 

8 

‘bgr-bfe1  ‘;A‘‘  ■<&»■  ‘tf  Ar  ‘rfAr  afe 


10 

- j  .£•  '-4.fr-'  '-4  b‘*rV,;,.^-  ■  ■  ,V  •  r,j\y  •,; ),r  '4.k 


■Jkps 


*  *  ?/  • 

•>'  5§* 

ro  ACS'  - 


>  r.  w  . 

l.  /J«<  fl.y.  cK 

..  ..  12 

■olcfe'  '-dife'  '-db*'  '-iiit'  '6h£'  'dAl'  'djfe 


life 


: 


v  '  *;  »y  •  •  t  - 

A  -  W 


rfr  ■  t:*: 

: 

14 


‘-6i£‘  ‘(jjK ‘flfe 


Figure  21.  The  first  16  individual  PIV  images:  position  22  (Dpt  998). 
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Figure  22.  Conditional  averaging  with  images  #2,  #3  and  #8:  position  22  (Dpt  998). 
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Vorticity,  BL  Pos  22,  Dpt  998 
rotated  with  images  #2,  #3,  #8 
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Figure  23.  The  vorticity  (y/Q)  and  cross  flow  (vz/QR)  averaged  with  images  #2,  #3  and  #8  for 
position  22  (Dpt  998)  before  and  after  the  rotation. 
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Figure  24.  A  comparison  of  the  centroid  and  convolution  of  X2  rotation  results  for  each  individual 
image,  #2,  #3  and  #8  and  its  averaged  image  for  position  22  (Dpt  998).  The  label,  -ea”  is  the 
conditionally  averaged  image. 
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a)  Vorticity  b)  Swirl  velocity 

Figure  25.  Cutting  axes  in  the  vorticity  map  and  the  associated  tangential  velocities  at  0°,  90°,  180°  and 
270°  for  position  17  (Dpt  1003). 
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Figure  26.  Cutting  axes  in  the  vorticity  map  and  the  associated  tangential  velocities  at  0°,  90°,  180° 
and  270°  for  position  22  (Dpt  998). 
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Figure  27.  Swirl  velocities  using  a  Vatistas  model  of  the  prototyped  images  for  positions  17  and  22. 
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